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OX) 1 Abstract 

In this paper, we develop a partner selection protocol for enhancing the network lifetime in 

o 

cooperative wireless networks. The case-study is the cooperative relayed transmission from fixed 
indoor nodes to a common outdoor access point. A stochastic bivariate model for the spatial 
distribution of the fading parameters that govern the link performance, namely the Rician K-factor 
and the path-loss, is proposed and validated by means of real channel measurements. The partner 



selection protocol is based on the real-time estimation of a function of these fading parameters, 
i.e., the coding gain. To reduce the complexity of the link quality assessment, a Bayesian approach 
is proposed that uses the site-specific bivariate model as a-priori information for the coding gain 



estimation. This link quality estimator allows network lifetime gains almost as if all K-factor values 

00 ' 

were known. Furthermore, it suits IEEE 802.15.4 compliant networks as it efficiently exploits 



the information acquired from the receiver signal strength indicator. Extensive numerical results 
highlight the trade-off between complexity, robustness to model mismatches and network lifetime 
performance. We show for instance that infrequent updates of the site-specific model through K- 
factor estimation over a subset of links are sufficient to at least double the network lifetime with 
respect to existing algorithms based on path loss information only. 
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I. Introduction 

Wireless sensor networks (WSN) are the enabling technology for home and building 
automation 0], O. The main obstacle in the development of WSNs is the cost of battery 
replacement, which becomes even more pronounced for indoor-to-outdoor (120) communica- 
tion with its larger transmit powers requirements. In this paper we aim at designing a medium 
access control (MAC) protocol for 120 WSNs that minimizes the maximum transmit energy 
so as to prolong the network lifetime 0. 

Transmit energy can be reduced by implementing advanced cooperative relaying strategies 
[IH, that efficiently exploit the inherent spatial diversity of a distributed radio channel. Se- 
lecting the partner for each node Q is the most crucial task in the coordination phase of 
any cooperative technique J6J. Partner selection is either based on instantaneous or average 
channel quality indicators of the links Q, (S). Additional knowledge on macroscopic features 
such as the network topology [0 and a parametric characterization of the fading channel, 
e.g., a path-loss model [10], can be also exploited. 

In this paper we focus on the transmission from several indoor static nodes in a single- 
floor office or factory to a common outdoor access point (AP), as outlined in Fig. Q] (top). 
Nodes are permitted to engage in cooperative transmissions by amplifying and forwarding 
(AF) [11 J the signals received from the partner nodes. We are interested to develop a partner 
selection scheme that maximizes the network lifetime under reliability and rate constraints 
that are identical for all nodes. We propose an original approach where the partner selection is 
aided by the knowledge of a site-specific multi-link stochastic channel model. The network 
employs a time division multiple access (TDMA) scheme inspired by the IEEE 802.15.4 
access protocol in "beacon mode" Ifl2l . as depicted in Fig. Q] (bottom). 

In selection cooperation based on instantaneous channel state information (CSI) has 
been introduced. For large networks such an approach requires the exchange of large amounts 
of CSI within the coherence time of the fading channel. To reduce the transmission overhead, 
here we are interested to employ long-term channel properties for the partner selection stage. 
This approach has been used by the pairing protocol in (S), which is based on the path loss 
values. This is a very practical choice, as e.g. in the IEEE 802.15.4 standard a received signal 
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strength indicator (RSSI) is available lfl3l . However, an indication of the randomness of the 
fading is also required for assessing the quality of a link. This randomness can be assessed 
by the Rician K-factor, as shown by physical fading channel modeling IT141 . 

Contributions of this work : (z) We design a MAC protocol where the AP assigns the 
partner and transmission resources to each node based on information about the path loss and 
the Rician K-factor, extending the method of JSJ. For this protocol we provide a performance 
assessment in terms of network lifetime, utilizing realistic I2I/I20 channel models. 

(ii) We propose an empirical, but analytically tractable, stochastic model for the char- 
acterization of the two fading parameters, namely, the path loss and the Rician K-factor 
for the links in an indoor network with fixed nodes. The two parameters are modeled jointly 
generalizing existing scalar models, e.g., lfl~4Tl. The so-called bivariate channel model is drawn 
and validated using the multi-link channel measurement data lfT511 . 

(Hi) For low-power operation of a cooperative WSN, we propose a procedure where (1) the 
K-factor is estimated on a small number of nodes with a slow update cycle: this distributed 
information is conveyed to the central coordinator in order to update the bivariate model 
parameters; (2) the estimation of the link-quality for each transmission is based on the 
local average RSSI and the regularly updated bivariate channel model as common a priori 
information. We compare the performance of this approach with the one where both the 
path-loss and the K-factor are permanently re-estimated for each transmission. 

II. System Model 

The scenario under study consists of N battery-powered indoor nodes that communicate 
with a common AP located outdoors. The nodes transmit during a communication session. 
In each session, the AP is acting as the centralized coordinator for assigning the cooperating 
partners, configuring the time-slot assignments and radio-frequency (RF) transmit powers. 
Complex base-band notation is used to model the wireless link between node i e {1, . . . , iV} 
and node j G {0, . . . , N}, with node j = referring to the AP. The received signal at node 
j is yij = hijXi + rij where hij is the frequency-flat complex channel coefficient, Xi is 
the symbol transmitted by node i with transmit power p i5 rij is additive symmetric complex 
white Gaussian noise with variance a 2 . The instantaneous signal-to-noise ratio (SNR) for 



August 11, 2011 



DRAFT 



3 



from node i to j is modeled as 

= (Pi/o*) \hij 2 ■ (1) 

The square fading envelope \hi t j\ is constant for the whole codeword duration (block fading) 
and varies from codeword to codeword according to the Rician distribution lfT6ll . such that 
Ljj = — (E [|^.jj| 2 ]) dB denotes the path loss in dB and K itj = (|E [/ijj]| 2 ) dB - (var [hij]) dB 
denotes the Rician K-factor in dB. The term path loss, used for indicating Ljj, includes the 
large-scale shadow fading and the static component of the small-scale fading as detailed in 
Sect. IIII-AI The block fading assumption motivates the use of outage probability and it is 
confirmed by channel measurements (see lfT5T0 . 

A. Medium Access Control and Node Coordination 

The transmission is organized into frames of duration Tp, further divided into N + 1 
subframes for time division multiple access (TDMA), as shown in Fig. Q] (bottom). A unique 
subframe of duration T s = T F /(N + 1) is assigned to each node by the AR The AP also 
provides the reference clock to all the nodes, the grouping decisions and access coordination 
(e.g., power and time slots allocation) through a periodic beacon transmission Ifl2l . 

Let be a pair of cooperating nodes^, to accommodate cooperative transmission each 
subframe assigned to any of these two nodes is further subdivided into two slots. As depicted 
in Fig. [0 for node i the first slot spans a fraction^ = /3 of the subframe duration and it is 
used to transmit the node i data. The second slot with duration (1 — /%) Ts = (1 — f}) Ts is 
reserved for helping the assigned partner node j. In a specular way, the first and second slots 
of the subframe assigned to node j span, respectively, the fractions 0j — 1 — ft (for node-j 
data) and (1 — fy) — (3 (for forwarding node-i data). The AP optimally combines the noisy 
replicas of the signals coming from the two nodes for data detection. 

'We assume that each node can cooperate at most with one partner. The extension of the analysis to grouping assignment 
with more than one partner is beyond the scope of the paper. 

2 To simplify the mathematical treatment, we assume that the slot partition ft can take any value in the interval < /3 < 1, 
although in practice this is constrained to a finite number of data units. 
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B. Outage Probability Modeling 

For point-to-point transmission the outage probability is P* t = Prfyy < 7$f]> where the 
SNR threshold 7^ = (2^ — l) /Y , R is the spectral efficiency measured in bit/channel use. 
The gap < Y < 1 can be varied by changing modulation/coding format and targeted bit 
error rate (BER) level (see [17] and references therein). The outage probability is assessed 
according to the models in lfT8ll . here adapted to the considered scenario. 

For the Rician fading model, the outage probability for node i communicating directly 
with any node j can be parametrized in terms of the so-called coding gain cy 

~ (2) 



= — T \ n i f\ ( Ty M ' ( 3 ) 



where m indicates that the equality holds asymptotically for high SNR|f|. The coding gain cy 
depends on the K-factor K it j and the path loss Ly according to: 

e{L itj )^d{K i}j )) 

For convenience of notation we introduce the function O(-) = 10^ '^ 10 with inverse = 
101og 10 (-) (recall that K i: j and Ly have been defined in dB). 

To model the performance of cooperative transmission of node i with the help of node 
j (here j ^ 0), we consider the AF relaying scheme as it has a simple architecture that 
facilitates practical implementation. Node j periodically overhears the signal transmitted by 
the partner node i and amplifies-and-forwards it towards the AP [fTTTl . The amplification is 
based on a variable gain approach, where the power amplification gain a, = pj/ [^('jij + 1)] 
is dynamically adjusted to the instantaneous SNR 7y. Note that the node must maintain a 
constant transmit power pj for the whole assigned subframe, i.e. also during the relaying 
phase, so as to avoid amplifier non-linearities by switching the power level. Recalling that 
the AP optimally combines the two noisy replicas of the source signal, the effective SNR 

/ 1 1 1 V 1 

7(«),o = 7i,o + — + + • (4) 

'Tightness of the approximation is verified for SNR large enough to guarantee sufficiently low outage probabilities 
(< 10 ')• 
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The outage probability for cooperative transmission 

P G C U T = Pr[7(u),o < 7 t 7 P ] - ~ ? - 7th „ , (5) 

with 7th° P = (2 i?//ft — l) /r. The spectral efficiency i? is now multiplied by l/fa to guarantee 
the same efficiency as for the non-cooperative case. Notice that this assumption implies that 
the low-power radio transceiver is designed to support multiple data rates lfT9l . Also, we select 
the same gap T for all rates. Terms P(ij),o, C(i,j),o and <%j) denote the effective power, coding 
gain, and diversity order for the cooperative link (i,j),0, respectively. For Rician fading, it 
is <%,j),o = 2, P(i,j),o = sfPifj , and 



c (m),o 



Cifi \Ci,j Cj,0 



(6) 



where the coding gains depend on the K-factors and the path loss values, as in ©. 



C. Energy Consumption Modeling 

The transmit power is designed so that the outage probability at the AP is lower or equal 
to p. The corresponding energy consumption for node i during one frame is derived below. 
No-cooperation - From ©, the transmit power at node i scales as 

P, « ^l 1 -. (7) 

Q,o p 

The average energy expenditure for node i can thus be modeled as 

Ei = Pl T s + E RX ^ + Ep, (8) 

J- F 

where is the energy consumption for receiving during the beacon slot only, and E P is 
the energy consumption for basic processing. 

AF cooperation - For a small enough outage probability p, it can be shown that the 
minimum transmit power levels for paired nodes are 



Pi = Pj = P(i,j),0 = /</3)cr /V 2 P, ( 9 ) 

with k((3) = max [(2^ - l) / (Tc {i>j)>0 ) , (2^^-^ - l) / (Tcy Afi )]. The subframe frac- 
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tion j3 for the z-th node message is selected so as to minimize P(ij),o : 




1 log 2 (A) 

2 8R 



(10) 



with A = \ + c '* o/c ^ . The notation « indicates that the scaling law is valid for log, (A) <C R 
(i.e., A ~ 1). The proofs of © and (flOl) are given in Appendix A. 



reserved for helping the partner j that experiences more severe fading conditions. On the 
other hand, if ^> max [c^o, c^o], the slots have equal length, j3 ~ 1/2. Notice that the 
choice (3 = 1/2 is relevant also because it models a practical system optimized for two data 
rates [fT9ll . i.e., the lower for no-cooperation and the higher for cooperation. 
The average consumed energy for node i cooperating with partner j is then 



with v AF > accounting for the energy consumption for the partner signal amplification. 



III. Wireless Channel Characterization 

In this section we present a framework for modeling and estimating channel quality metrics 
for the links of the cooperative network. The final aim is to provide metrics to be used at 
the MAC layer for relay selection. 

In Sect. IIII- A1 we propose a stochastic general model of the two fading parameters (L^j, K^j), 
which includes as particular cases some models previously proposed in the literature (see 
e-g- El, EOlO . In Sect. IIII-Bl the model is validated and discussed using experimental data 
collected by an I2I/I20 measurement campaign, with indoor nodes deployed over an office 
environment. The study will also provide a reliable simulation environment for assessing the 
performance of partner selection algorithms in Sect. |V] Given the a-priori knowledge of the 
model, in Sect. IIII-CI we propose a Bayesian method for the estimation of the coding gain 
Cij starting from the path loss observation. 



Notice that using (flOl) for A > 1 (i.e., c i > Cjo), then < 1/2 as the largest slot is 



= P(iA0 T s + E RX (2 - /3)-i + (l + v AF )E P , 



(11) 



The transmit power P(i t j),o for source and relay nodes are chosen as in ©. 
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A. Bivariate Model for the Large-Scale Fading Parameters 

We consider a fading channel hij between any two nodes of the static network. Link 
indexes will be omitted when not needed, to simplify the notation. It is important to 
highlight that, differently from mobile scenarios, here the deterministic channel gain fx h = 
E [hij] accounts for the effects of fixed scattering/absorbing objects which determine the 
so-called static multipath component. On the other hand, temporal fading - with variance 
a\ = var [h i: j] - is due only to some moving scatterers/absorbers in the environment. As 
assumed in Sect, [nl this results in a Rician distribution with parameters L = — 9~ 1 (af l + \fih\ 2 ) 
and K = O^d/i^ 2 /al) (recall that #(■) is the inverse of the transformation to dB). The 
multipath configuration changes rapidly with the node locations, thus leading to fast variations 
of the Rician factor K and the path-loss L over the space. The objective of this section is 
the definition of a model to describe the variations of such parameters from link to link. 

Let x = [i^ L] T be the vector collecting the two fading parameters for a generic link. 
We model the variations of x according to a bivariate Gaussian random variable (bivariate 
model). We assume that x is Gaussian distributed, x ~ A/"(/x x , C), with a mean /jl x (D) = 
[/ik(-D) /xl(-D)] t depending on the link distance D and a spatially invariant covariance matrix 



C 



<7 K ip<J K (JL 
^2 



(12) 



(pa K a L 

The correlation coefficient tp = E [{K — /xk(-D)) (L — /(&}<.<? l) models the mutual 

dependence between the K-factor and the path loss experienced over the same link. A metric, 
that will be relevant in partner selection analysis in Sect. [IV] is the probability density function 
(pdf) of the K-factor conditioned on the path loss, p (K\L) = Af(/J>K\u °k|l)' w ^ m mean 

fi K]L = E [K\L] = fi K (D) + — (L — fi L (D)) (13) 

and variance a^ L = var [K\L] = (1 — tp 2 ) a K- 

The parameters of the bivariate model can be tuned for different propagation scenarios, 
e.g. 121 and 120, as done in the following. 
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B. Experimental Calibration of the Bivariate Model and Analysis of Spatial Coherence 

In this section, we describe the calibration of the model parameters /U X (D) and C on the 
121 and 120 propagation scenarios in Fig. [TJ We use the multi-link channel measurements 
lfT51l at 2.45 GHz, but also existing models whenever the data is not sufficient. It is important 
to mention that the fading in the experiment is caused not only by walking people but also 
by moving metallic objects. This contributes to determine harsh 120 propagation conditions. 
The 70 MHz band of the channel is divided into 60 subbands, each corresponding to a flat 
fading subchannel. As mentioned in Sect. IIII-AL in contrast to mobile scenarios, the small- 
scale temporal fading has a different origin compared to the small-scale fading in the spatial 
and spectral domains. Hence, K^j and L^j are estimated independently in each subband. 

121 channel model - The vector function /i x (-D) is estimated by performing linear least 
squares regressions of K and L over the corresponding distances D in logarithmic scale. 
The covariance C is then obtained by computing variances and covariances on the sets 
of data {K — fx K (D)} and {L — /i L (£>)}. The norm of the error between the theoretic 
cumulative density function (cdf) of {K — fi K (D)} and {L — fi L (D)} and the empirical one 
is 0.04, showing a very good fit. This is apparent in Fig. O where the equidensity contour 
lines of the bivariate distribution of (K — jj K (D), L — fx L (D)) are shown together with the 
respective measured values. 

120 channel model - According to the urban micro-cell scenario B4 in ETTl . the propagation 
on the link (i, 0) is modeled as the combination of three main contributions: (0 the indoor 
propagation from the node to the nearest wall to the BS (i, Wall); (if) the propagation through 
the wall; (Hi) the outdoor propagation from the wall to the BS (Wall, 0). The overall link 
path loss Lj o and K-factor K i are modeled as 

Lifl = £j,Wall + -^Wall + Avall.O) K i0 = -K^Wall + -^Wall.C)- (14) 

Notice that the wall contribution has no effects on K^. The value chosen for the wall 
contribution is L Wall = 14dB (neglecting the angle of the propagation path with respect to 
the wall [21]). The specific value will not affect the performance comparison of the algorithms 
presented in Sect. |V] The outdoor parameters (Lwaii,o, -Kwaii.o) modeling adheres to ll2~0l (4) 
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and (9)]. The indoor parameter L^waii is modeled according to ETl Table 4-4]. The bivariate 
model wall j ^Q,Waii) an d the corresponding indoor correlation value are obtained by a 
simple manipulation of the model for L^waii in [21, Table 4-4], using the line slope and the 
error variance of the least squares linear regression of K over the corresponding L in the 
available 120 measurements. Numerical details are provided in Table 1. 

The normal distribution of K \ L has been assumed by several studies in the literature 
(see, e.g., Il20ll . Iffifl . |[T4l ) to model the residuals of the linear least square regression 
of K over L. Here, instead, we have directly tested the bivariate normal distribution of 
{K — hk(D)} and {L — ^(D)} via a 2-D fitting of the cdf. In this way, we have also 
highlighted the strong (negative) correlation between {K — /x K (L>)} and {L — /x L (D)}, also 
observed experimentally in our study (see the values detailed in Table 1). In contrast to 
line-of-sight mobile scenarios, in our case the static multipath component in jj, h mainly 
contributes to the negative correlation ip. Consider e.g. two links and (i, j 2 ) ? with 

nodes ji and j 2 closely located, i.e., with distance in the order of the carrier wavelength A 
(Di w Dij 2 ). Even if the two links are likely to experience the same channel shadowing, 
average fading conditions // X (D), and degree of temporal variations^] er^, the two links still 
experience different multipath configurations and thus different values for the static multipath 
component in fx h . The variation on /i h affects with opposite sign the path-loss and the K- 
factor, thus, a strong negative correlation is observed. 

Spatial coherence of the models - We use the available measurement data to empirically 
assess the spatial coherence of the 121 bivariate model. We first estimate the model from the 
measurements over 16 links during the first 1/3 of the total experiment duration (32s). The 
model is again estimated from a disjoint set of 16 links measured during the last 1/3 of the 
total experiment duration. We observe that the model based on the first subset of links can 
predict with extremely high accuracy the mean values /x x (D) as modeled from the second 
subset of links. Fig. [3] shows the measured (L — fi L (D) , K — fi K (D)) values for the two 
subsets of links and the corresponding bivariate Gaussian contour lines. The depicted bivariate 
models are strongly matching as the covariance matrices C are practically the same. Due to 

4 This has been observed in the frequency domain also in 1 14 1. Notice that the small-scale spatial and frequency-selective 
fading are caused by the same mechanism. 
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the common underlying physical mechanism, it is realistic to assume that similar results are 
valid for the indoor propagation model (L^waib ^Q.Waii) m me 120 case. Finally, the outdoor 
(£waii,o, -Kwaii.o) values, similarly to the wall penetration loss L Wa n, are not changing over the 
links and time in the considered 120 scenario (so as to model a common stationary outdoor 
propagation scenario). As analyzed in Sect. IV-Cl these conclusions suggest that the model 
parameters can be extracted efficiently employing just a small subset of nodes. 

C. Link Quality Estimators 

In indoor environments, it is likely to incur into a link that exhibits both a large average 
RSSI and, yet, a high packet loss rate lfT3ll . In these scenarios, accurate link quality estimation 
should include a measure of the fluctuations of the received power. Based on the physical 
fading channel modeling, several works propose ways to assess the randomness of a link via 
the online estimation of its K-factor. However, it is not trivial to predict the physical layer 
performance, or even to draw MAC layer decisions directly using the estimated K-factor. 

Let us instead consider the problem of the estimation of the coding gain ©. It is convenient 
to re- write the coding gain in dB scale as: 

c = T9(K) -9~ l [1 + 00*0] ~L, (15) 

V V ' 

where T = 101og 10 (e). The term s(K) measures the additional information provided by the 
coding gain in Rician fading (9(K) > 0) with respect to path-loss information only. 

Most commercial RF transceivers designed for low-power wireless applications Il22l . ffT9ll 
provide information about the received signal strength from which the path-loss L can be 
easily inferred [fT3l . A straightforward method to estimate c would then be to calculate also 
the K-factor and to adopt the formula ©. We refer to this estimate as direct estimate cfj. 
Accurate, but complex, estimators are proposed, e.g., in ll23Tl . The most accurate estimator 
requires the knowledge of the fading phase. On the other hand, for non-coherent transmissions 
as in current implementation of IEEE 802. 15.4 radios, the ratio between the squared mean and 
the variance of the fading envelopes can be also used to estimate the K-factor. Nevertheless, 
the estimation of small K-factor values (K < 3dB) becomes very inaccurate [1241 . 

Here, to reduce the estimation complexity, we propose to estimate c only from the path 
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loss L by exploiting the a-priori information on the statistics of c\L derived from the bivariate 
model. We assume the stationarity and the perfect knowledge of the parameters /i x (-D) and 
C introduced above. It is important to stress that the model knowledge is available only if the 
nodes are cooperating to exchange the information required to extract the model parameters f]. 
The assumption of model knowledge is therefore realistic in the considered network. 

Bayesian estimation of the link quality c based on the observation of L requires the 
computation of the a-posteriori pdf p(c|L). This pdf can be approximated by observing that 
for T6(K) > 6- 1 [1 + 9(K)] it is q(K) w T6(K) and the link-quality indicator <Q3) reduces 
to c ~ T9(K) — L. Recalling that K conditioned on the observed L is Gaussian distributed, 
K ~ JV ^A*k|l ; ^kilJ' it follows mat c is shifted log-normal distributed with pdf ll25Tl : 

p( c |L) ~ ^27r4 |L ] _1 exp (-^- [6' 1 (£±£) -/i K | L ] 2 ) u(c + L), (16) 

where u(x) is the unitary step function. The approximation is corroborated by a numerical 
analysis in Fig. HI where the true and the approximated pdfs are shown for two values of the 
path loss L of one 121 link, simulated according to the model calibration in Sect. IIII-BI 

Based on the a-posteriori pdf, we derive the maximum a posteriori (MAP) estimator of the 
link quality c. The MAP estimator, c MAP = argmax [p(c|L)], is approximated by the mode 

c 

value of the shifted log-normal distribution (fT6l) which yields Il25ll : 

c MAP « T9 ^ K | L -^j-L. (17) 

A minimum mean square error (MMSE) estimator c MMSE is also provided in Appendix 
B. The MAP and MMSE estimators are depicted, by empty markers for L = 63dB, and by 
filled ones for L = 58dB in Fig. |H together with the corresponding pdfs. The figure shows 
that the approximations of the pdfs and of the estimators are tight to the exact ones, and 
behave similarly with varying path loss value. 

IV. Partner Selection Strategies 

In this section, the problem we tackle is how to select the partner for the nodes that 
are tasked to communicate to the outdoor AP. The final goal is to minimize the maximum 

5 The distributed estimation of the 121 and 120 models is beyond the scope of the present work. 
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energy consumed among them. The AP selects the partners (relays) based on long-term link 
quality metrics and not on the instantaneous channel gains, as proposed and studied in the 
literature (e.g., in Q). The optimal min-max pairing is found in Sect. IIV-AI that allows also 
for configurations with un-paired nodes. To lower the complexity and the amount of signaling, 
a worst-link-first (WLF) algorithm is then described in Sect. IIV-BI Although the structure 
of the WLF algorithm is simple and resembles the one in [[8]| - with a modification for odd 
number of nodes - the proposed performance analysis diverges substantially for the I2I/I20 
fixed links considered in our experimental scenario (see Sect. |nl]). Given that the distributed 
wireless links can be modeled by Rician fading with different K-factors, the first key idea 
is that the AP uses the coding gain and not the path loss as decision metric in the WLF 
algorithm. Secondly, the Bayesian estimators derived in Sect. IIII-CI provide a representation 
of the link quality to be used by the AP to finalize the partner selection. 

A. Problem Definition and Optimal Solver 

We define the set of candidate pairing sets V, such that one set £ E V contains up to 
L-/V/2J disjoint pairs of cooperative nodes: £ = (k, h), (/,<?)} ■ All the non-paired 

nodes belong to the set of single nodes S$ = {q, s, . . . , z}, such that 2 |£| + |<Sf| = N 
(where |-| denotes the cardinality of the set). Given the candidate pairing set £ and the 
corresponding single node set S%, the maximum energy consumed by a node in the network is 
-E max (£) = max[maxE^ , max£ ? ], where E™* fi = max[E^j) j0 , i?(j,i),o] is the maximum 
energy for the pair for a given relaying protocol, with E^j^ defined by (TTTT) . The 
optimal pairing £ is the solution to 

| = argmin£ max (0 , (18) 

£,ev 

where we assume that all nodes have the same rate R and outage probability p constraints. 

The problem (fl"8l) can be formulated as a special case of the weighted matching prob- 
lem on the non-bipartite graph Q = (X,£) [26 J . The set of vertices X corresponds to the 
set of nodes {1, . . . , N}, which are fully connected by the set of undirected edges £ = 
{eij : E X) & (i < j)}. The loops ej J= , can be regarded as edges e i i, where the virtual 
vertex i of the extended graph is connected only to i. The weights w (e^i) = E™* and 
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w (eij = i) = Ei are associated to all the edges and loops, respectively. The optimal pairing 
algorithm removes at each iteration the maximum weighted edge of the extended graph (as 
done in [271 for a max-min problem on the bipartite graph) and checks the existence of a 
weighed matching solution in the remaining graph using Gabow's algorithmic [26, Ch. 11], 
which was instead proposed in [[81 for minimizing the sum of the energies consumed at the 
nodes in one iteration only. It can be shown that the above algorithm reaches the solution in 
O (N 5 ) computational time. Notice that the algorithm is centralized and requires the AP to 
know all the inter-node link qualities for computing cuj\ . 

B, Worst-Link-First Coding-Gain Based (WLF-CG) algorithm 

The WLF algorithm is a suboptimal protocol for node pairing that allows to reduce the 
complexity to O (N 2 ). The conventional WLF algorithm (referred to as WLF path-loss-based, 
WLF-PL) is based on the information of second order statistics of the fading link [8|, i.e. 
the path loss L iy j. We propose a WLF method based on the coding gain (WLF-CG). 

Link quality estimation - Before pairing decisions can take place, each node i locally 
acquires an estimate of the link qualities for all 121 links Cjj to the candidate partners and 
the 120 link towards the AP Cj >0 . Link quality estimator options cfj, c^ AP , and c^mse (in 
dB) have been discussed in Sect. IIII-Cl The WLF-PL algorithm uses dij = —L^. 

Protocol structure - The algorithm is composed of two phases: 

1 ) Candidate partner set discovery - Each node % performs link quality estimation c ij0 for 
the 120 channel exploiting the periodic transmission of beacon subframes from the AP as 
probing signals to be used for channel parameter estimation. The link qualities Qj measured 
from all neighboring nodes over 121 links are estimated based on the signals overheard from 
the potential partners. The difference qj — Cjo is compared at each node i to a common 
threshold r in order to guarantee the condition 0(ci t j) S> 9(ci t o). If — c^o > t then node 
j becomes a candidate partner for node i. Given the 120 link quality estimation c i)0 , the 
threshold r is centrally designed such that the probability of finding no candidate partners 
among the N — 1 potential candidates [X^Li^ [1 ~~ P r ( c m > ^i,o + r l-^ij)] is small enough, 
where Pr (qj > c i + T \Lij) can be estimated using (TT6l) . The evaluation of the optimal 

6 The Hungarian method, which is tailored for the bipartite graphs, is instead considered in 1271 . 
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value for the threshold r is carried out in the case-study outlined in Sect. [V] The candidate 
partners set Cp(z) = {j : cy — c i0 > r} is finally communicated to the AP from each node 
i, using the assigned subframe. 

2) Assignment algorithm at the AP - At each iteration the AP selects the worst-uplink 
node i with link quality c i;0 < Ck,o >V£; ^ i, and, if possible, assigns it to the best-uplink 
candidate partner j, such that 

j = arg max c,-n (19) 

Nodes i and j are paired and disregarded in the next iterations, unless C-p{%) is empty. In the 
latter case, node i is left un-paired in the final configuration. 

For an odd number of nodes N the preliminary step is to leave the best-uplink node 
un-paired, then the partner assignment follows as above for the remaining N — 1 nodes. 

V. Experimental Assessment of Partner Selection in the 120 Environment 

In what follows, we provide numerical simulations on the performance of the partner 
selection algorithms (Sect. HVb and of the Bayesian link quality estimation methods (Sect. 
IIII-CI) . The network setup is outlined in Fig. [H N nodes are randomly distributed in a 
25m x 25m indoor environment, while the AP is placed outdoors 50m away from the 
nearest wall. For each random topology of the network, K-factor and path loss values 
are generated independently for all links according to the respective 121 and 120 models 
detailed in Sect. IIII-BI The performance results are expressed in terms of lifetime gain 
E[^ max (£ ref )]/E[^ max (£)] ; defined as the ratio between the maximum energies (averaged 
over 5 x 10 4 random topologies) consumed among the nodes according to two different 
pairing strategies: the reference protocol with pairing solution £ rc f and the proposed protocol 
with pairing solution £. Note that the largest values of E max (^ re{ ) and of i? max (£) highlight the 
lifetime gain in scenarios where the 120 channel conditions are worst. In the examples, the 
energies consumed for reception E^x. and for basic processing E v are neglected. According 
to the 120 modeling in Sect, [nil it is likely that cy max [c$ , Cj,o\- From (flOl) . the optimal 
choice for slot duration is (3i = /3j = 1/2. We set T — 1, as this assumption has no relevant 
impact on the performance comparisons in the examples. The target outage probability is 
p = 10~ 3 for all the nodes with spectral efficiency R = 1 bps/Hz. 
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In Sect. IV-A[ we provide numerical simulations for the WLF-CG lifetime performance. 
In Sect. IV-B[ the impact of the proposed link quality estimation is evaluated and compared 
to the case where a noisy estimation of the K-factor is obtained from training. Finally, in 
Sect. IV-Cl the proposed pairing algorithm performance is discussed in presence of imperfect 
modeling of the channel. 

A. WLF-CG Protocol Performance 

We first evaluate the performance of the WLF-CG algorithm assuming that the link qualities 
Cjj are perfectly known at the respective nodes i — {1, . . . , iV} and j — {0, . . . , N}. Fig. [5] 
shows the lifetime gain of AF cooperation compared to no-cooperation, i.e. £ ref = is the 
empty pairing set and £ is the pairing set obtained according to a partner selection strategy. 
A random pairing strategy is also considered for comparison where all nodes are disjointly 
paired with a random choice of the partner. For odd N, the optimal pairing algorithm in Sect. 
IIV-AI the random pairing strategy, the WLF-PL algorithm, and the proposed WLF-CG are 
compared. The candidate partner conditions (see Sect. llV-Bl) cy c i;0 and Ly <C L i for the 
WLF-CG and the WLF-PL, respectively, are almost always guaranteed for each pair of nodes, 
in particular it is Pr (cy > q i0 + t) ~ 1 and Pr (Ly < L^o — r) ~ 1 for r < 30dB. For 
the propagation environment under consideration, the simulations show that partner selection 
performance get worse when choosing r > 40dB. The exploitation of the knowledge of the 
K-factor is revealed crucial: the WLF-CG algorithm increases the lifetime from a factor of 
20 for iV = 3 to a factor of 2 for iV = 55 compared to the WLF-PL. This results from 
the fact that the WLF-CG algorithm allows for a more efficient exploitation of the available 
diversity, as if the optimal algorithm in Sect. IIV-AI were applied. The remarkable gains over 
no-cooperation and over the random pairing denote a large degree of spatial redundancy 
provided by the multi-link channel, i.e., the path loss and K-factor values exhibit significant 
variations over the space. 

B. Impact of Link Quality Estimation on WLF-CG 

The optimality of the WLF-CG pairing strategy relies on the accuracy of coding gain (link 
quality) estimation. This motivates a closer study to identify the most suitable estimator and 
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to quantify the benefits provided by the knowledge of the distributed channel model. Here, 
we assume that only the path loss L is perfectly known. 

Link quality estimation from the path loss cy = {c^ AP , c^ MSE } - The WLF-CG algorithm 
can capitalize from the available site-specific channel characterization. Fig. [6] shows the energy 
gains of cooperation over no-cooperation for varying number of cooperating nodes N. The 
estimators c^ AP and c^ MSE are used, exhibiting equivalent performance: notably, the lifetime 
gain over the WLF-PL ranges between factors 2 (N = 15) and 20 (N = 3). These gains are 
similar to those obtained in the case where cy are perfectly known (see Fig. |5]). Thus, the 
proposed Bayesian estimation of c is revealed useful to guarantee significant lifetime benefits 
in the considered network settings. 

In Fig. |7] we evaluate how the WLF-CG lifetime gain scales with the correlation tp between 
the path loss and the K-factor for N = 10. The correlation ip (assumed equal for both the 121 
links and the indoor component of the 120 links) varies arbitrarily from —1 to 0, whereas 
the other channel parameters conform to the values in Table 1. When the link quality c is 
known, the lifetime gain increases as cp gets closer to 0: this is the case where the knowledge 
of the K-factor becomes in theory most beneficial. Instead, the WLF-CG algorithm based on 
c^ AP , and c^ MSE is shown not to be sensitive within realistic variations —0.8 < cp < —0.5, 
as it improves always by a factor 2.5 the lifetime obtained by employing the WLF-PL. 
It is interesting to observe that Cy AP performs better than c^ MSE for ip > —0.4, as it is 
more conservative in predicting the worst-uplink link quality, which dominates the lifetime 
performance. Indeed, it can be verified that c^ AP < c^ MSE . 

Link quality estimation from both path loss and K-factor dij = Cy- - The estimator of 
the K-factor can be obtained from the complex fading realizations or the corresponding 
squared envelope values (acquired from RSSI measurements) ll23Tl . As discussed in Sect. 
IIII-Cl depending on the choice of the estimator, various trade-offs between accuracy and 
complexity can be obtained. Here, we prefer not to consider any specific estimator K, but 
we rather model the estimation noise 9 ^A'j — 9 (K) as zero-mean Gaussian, where the 
pdf is truncated in order to keep 9 (^K^j positive. Fig. [8] shows the lifetime gain of the 
WLF-CG algorithm with the direct estimate Cjj = cf^ over the WLF-PL, with varying root 
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mean squared error (RMSE) a AK = wE 



. Remarkably, for N = 10 and 



efk^j-e (k) 

o'ak < 5dB the lifetime is at least doubled. Notice that c^ AP outperforms c^ MSE , resulting 
in lifetime gains as if Cy were employed with ctaa' = OdB, and 3dB, respectively. The 
estimator of cfj is revealed robust for partner selection for a small number of cooperating 
nodes N. Instead, for larger number of nodes iV > 30, the lifetime is doubled only for 
a isk < — 10dB, while the WLF-CG is even outperformed by the WLF-PL for a& K > 6dB. 

C. Impact of Imperfect Model Knowledge 

We consider a protocol where a set of Ni nodes is used in a prior communication session to 
update the channel model, as shown in Fig. \T\ In a later communication session, N 2 different 
nodes are then tasked to transmit to the AP. The impairment of the model knowledge is due 
to the limited number of the (K, L) regression points used in the prior session to estimate 
the model, i.e., (iVf — N\)/2 points for the 121 model and N\ points for the 120 model. Thus, 
the accuracy of the proposed Bayesian link quality estimation improves for increasing N\, at 
the expense of spectral and computational power efficiency. This trade-off is assessed in the 
following. We focus on the MAP estimator c^ AP , that was shown to have better performance 
compared to c^ MSE in the above evaluations. 

Fig. [9] shows the pdfs of the 120 coding gain estimation absolute error in dB A c = |cj )0 — c i>0 \ 
for the estimator c^ AP with perfect (solid lines) and imperfect model knowledge derived 
via A r 1 = 7 nodes (dashed lines) for K ifi = {OdB, 7.8dB}. Notably, the absolute error 
pdfs are very similar, although the bivariate model is estimated by only 7 spatially distinct 
measurements of (K^o, L it0 ). As expected, the mean points of A c in the considered range of 
K-factor values, i.e., >0dB and <7.8dB, are smaller than the absolute error of the link quality 
estimator used in the WLF-PL, i.e., ${K) = c^o + L ii0 (depicted by the cross markers). 

Fig. [TO] plots the network lifetime gains of the WLF-CG that uses Qj = cf]^ F over the 
WLF-PL with varying number of the nodes Ni employed for the parameter extraction. The 
performance with perfect model knowledge is also depicted as upper bound. It is revealed 
that, for 3 < N 2 < 11, the proposed WLF-CG algorithm outperforms the WLF-PL, by at 
least doubling the network lifetime for N\ > N 2 (as highlighted by circle markers). 



August 11, 2011 



DRAFT 



18 



VI. Concluding Remarks 



We have designed an efficient partner selection protocol for cooperative wireless access 
from fixed indoor nodes to an outdoor AP through AF relaying. Given that the considered 
links can be modeled by Rician fading, we have proposed a partner selection algorithm that 
adopts the coding gain - a function of the path loss and K-factor - as link quality metric for 
pairing the nodes. The proposed approach implies an additional computational cost due to 
the estimation of the K-factor, but provides a network lifetime increase by factors ranging 
from 2 to 20 compared to the conventional algorithm based on path loss information only. 

Analyzing measurement data from a channel sounding campaign at 2.4GHz we were able 
to characterize the path loss and the K-factor with a Gaussian bivariate model. From this 
bivariate model a novel link quality indicator is derived that does not require a permanent re- 
estimation of the K-factors. This is the Bayesian estimate of the coding gain, where the path 
loss is the observed variable (inferred through RSSI readings) and the channel model is the 
a-priori information. The novel metric improves remarkably the partner selection performance 
almost as if full knowledge of the K-factors were available. 

Furthermore, the analysis on the measurement data reveals the high degree of spatial 
coherence of the model. Hence, we have proposed a protocol with two phases: (i) for a long 
time period (within the stationarity time interval of the model) low-complexity communication 
sessions take place, where the link qualities are inferred through path loss measurements and 
the a-priori information of the channel model; (//) for a short time period, a more complex 
communication session occurs, where a set of nodes estimate also the K-factors in order 
to update the site-specific channel model. Numerical results show a good trade-off between 
performance and robustness. In particular, the proposed protocol allows to double the network 
lifetime compared to the conventional algorithm, also in presence of modeling mismatches. 



For Rician fading where P(ij),o = y/PiPj and d(i,j),o = 2, the target outage probability p 
constrains the powers pi and pj over the two slots such that repetition based coding prescribes 
that pi = p and 1 - (3j = (3: 



Appendix A 




c r 2 (2«/( 1 -«-l) 



rc U,i),0 



{PPi) 1 , Pi>\ 




(PPj) 



-1 



(20) 
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Recall that the gap T can be designed to be the same for the communication both from i and 
from j. By minimizing the maximum over p^ and pj, the simple power balancing solution p^ = 

p 3 = K0)<j 2 /y/2pis found where k(/3) = max [(2 R ^ - l) / {Tc (iJ)fi ) , (2 R ^ 1 -^ - l) / (r C(j - 
and /3 = argrninfi;(/3) is solution to ^2 R ^ — cy^o = \2 R ^ 1 ~^) — C(jj) j0 , therefore 

2^-l = Af2 fi /M)-l) (21) 



l+Cj,o/cj,j 

such that for log 2 (A) <C -R it is 

R R 



where A = w c ''°, . The solution to (12TI) is now approximated for large enough rate R 



~log 2 A. (22) 



/9 1-/3 



Now by letting /3 = | — £ with £ small enough, the solution to (1221) is £> 1o | 2 ^ a ^ . 

Appendix B 

The minimum mean square error (MMSE) estimator c MMSE can be approximated as 

a MMSE = E[c\L}^E[T9(K)\L\-E[m8ix(K,0)\L}- L = 

- ™ (** + #) - «*« (-») " ^ «P (^) - 4 (23, 

where Q (•) is the Q-function. In (1231) we use the approximation <?(-rV) — T9(K) — max (if, 0) 
ViJll where E [max (K, 0) |L] = f Q °° ifp [AT|L] diT. 
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Indoor-to-outdoor networking: experimental scenario 

Office (Stanford 2008 measurements [15]) 
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to derive a site-specific channel model 
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Figure 1. Top: general indoor networking scenario and the specific 120 office radio measurement plan in 1151 : indoor 
nodes are allowed to support the estimation of a site-specific stochastic channel model (1), and to engage in cooperative 
transmission to the access point (2). Bottom: TDMA framing structure inspired to the IEEE 802.15.4 beacon mode. 




Figure 2. Measured 121 path loss and K-factor values minus the respective distance-dependent means. The equidensity 
contours of the zero-mean bivariate Gaussian distribution are depicted at one and two standard deviations from the mean, 
containing 62% and 98% of the points respectively. 



Table I 

Parameters for the Bivariate Model (fT2l ) - see Sect. IIH-BI 





Bivariate Model (1121) Parameters 




f^KpU) = 16.90 - 10a K log 10 p| m ) 
^ X = 40.4 + 10a L log 10 (D\ m ) 
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(i,j)^0; D = Dij 


[«k = 0.53, = 1.75] 
|V9 = -0.66 
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^ = 135.78 + 10« L log 10 PU) 

[ok = 0.45, a L = 3.89] 

c = -0.25 

1 o-k = 7.5dB, a h = 7.9dB 


Indoor- to-Outdoor (120), 
Indoor propagation if iiWall and L iiWall : 
(i ^ 0, wall); D = A.waii 
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= -0.74 
1 o-k = 5.7dB, cr L = 7dB 
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Figure 3. Measured 121 path loss and K-factor values minus the respective distance-dependent means for the first and for 
the second subset of links, estimated at two different time after 10 seconds, respectively. 
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Figure 4. Conditional pdf p(c\L), both exact (simulated) and approximated (analytical). The respective link quality 
estimators c^/" p and cJj VISE are also depicted for two values of observed path loss L (empty markers for L = 63dB and 
filled markers for L = 58dB). 
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Figure 5. Cooperative transmission lifetime gain over no-cooperation with different partner selection strategies and varying 
number of transmitting nodes N (odd values). The WLF-CG uses threshold values r = {30dB, 40dB} for Cjj — Cj,o in 
the candidate partner set discovery phase. The conservative choice r = — oo is used for the WLF-PL. 
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Figure 6. The lifetime gains of the WLF-CG algorithm based on c^ AP and c^ MSE over no-cooperation, compared to 
that of WLF-PL, with varying number of transmitting nodes N (odd values). 
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Figure 7. Lifetime performance of the WLF-CG algorithm, also with MAP and MMSE estimators, compared in percentage 
to the one of the WLF-PL (100% means equal performance). The correlation ip between path loss and K-factor as defined 
in d!2t is varying and is assumed equal for both the 121 links and the indoor component of the 120 links. As practical 
reference, the cross markers highlight the values <p = —0.66 and <p = —0.74, i.e., the correlation observed in the 121 and 
in the indoor component of the 120 measurements, respectively. 
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Figure 8. WLF-CG lifetime compared in percentage to WLF-PL, with varying K-factor estimation MSE. The lifetime 
performance of the WLF-CG algorithm based on c^ AP and c^j MSE are also marked on the respective curves for comparison. 
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Figure 9. 120 coding gain absolute estimation error |Ac| = |cj,o — Cj,o| pdf for the estimator c^o = P with perfect 
model knowledge (solid lines) compared to that with the model estimated via Ni — 7 points (Z/i,o,Jfi,o) observed in the 
previous communication session (dashed lines). Also the estimation error s(K) — a t o + Lj,o in the WLF-PL is depicted 
for the considered values of the K-factor K it a = {OdB, 7.8dB}, i.e., 9(Ki t o) = {1, 6}. 
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Figure 10. WLF-CG lifetime with MAP link-quality estimator cjj compared in percentage to WLF-PL, with varying 
number of nodes Ni involved in the the estimation the bivariate channel model during the previous communication session. 
The lifetime gains of the WLF-CG algorithm based on c^, AP with perfect model knowledge are also depicted with dashed 
lines as upper bound. Performance for different values of communicating nodes N2 in the current session are depicted. The 
circle markers highlight the performance for Ni = N2. 



